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Abstract. Let X be a compact manifold with boundary, and g a scattering 
metric on X, which may be either of short range or 'gravitational' long range 
type. Thus, g gives X the geometric structure of a complete manifold with 
an asymptotically conic end. Let H be an operator of the form H = A + 
P, where A is the Laplacian with respect to g and P is a self-adjoint first 
order scattering differential operator with coefficients vanishing at dX and 
satisfying a 'gravitational' condition. We define a symbol calculus for Legendre 
distributions on manifolds with codimension two corners and use it to give a 
direct construction of the resolvent kernel of H, R{a + iO), for a on the positive 
real axis. In this approach, we do not use the limiting absorption principle 
at any stage; instead we construct a parametrix which solves the resolvent 
equation up to a compact error term and then use Fredholm theory to remove 
the error term. 



1. Introduction 

Scattering metrics are a class of Riemannian metrics which describe manifolds 
which are geometrically complete, and asymptotically conic at infinity. We consider 
manifolds which have one end which is diffeomorphic to y x [1, oo)r, where K is a 
closed manifold, and is metrically asymptotic to dr^ + r'^h, where /i is a Riemannian 
metric on y, as r — > oo. The precise definition is given in Definition ^.1| below. 
Examples include the standard metric and the Schwartzschild metric on Euclidean 
space. 

In this paper we give a direct construction of the outgoing resolvent kernel, 
i?(cr + iO) = {H — (cr + iQ))~^ , for a on the real axis, where if is a perturbation of 
the Laplacian with respect to a scattering metric. The incoming resolvent kernel, 
i?(cr — zO), may be obtained by taking the formal adjoint kernel. 

The strategy of the proof is to compactify the space to a compact manifold X 
and use the scattering calculus of Melrose, as well as the calculus of Legendre distri- 
butions of Melrose-Zworski, extended by us in . The oscillatory behaviour of the 
resolvent kernel is analyzed in terms of the 'scattering wavefront set' at the bound- 
ary. Using propagation of singularity theorems for the scattering wavefront set leads 
to an ansatz for the structure of the resolvent kernel as a sum of a pseudodiffer- 
ential term and Legendre distributions of various types. The calculus of Legendre 
distributions allows us to construct a rather precise parametrix for the resolvent in 
this class, with an error term E which is compact. Using the parametrix, we show 
that one can make a finite rank correction to the parametrix which makes Id +E 
invertible, and thus can correct the parametrix to the exact resolvent. 

As compared to the method of 0|, where the authors previously constructed 
the resolvent, the construction is direct in two senses. First, we write down rather 
explicitly a parametrix for R[a + iO) and then solve away the error using Fredholm 
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theory. In Q, by contrast, the resolvent was constructed via the spectral measure, 
which itself was constructed from the Poisson operator. Second, we make no use of 
the limiting absorption principle; that is, we work directly at the real axis in the 
spectral variable rather than taking a limit R{a + ie) as e 0. We then prove a 
posteriori that the operator constructed is equal to this limit. 

Let us briefly describe the main result here. We consider an operator H of 
the form H — A + P acting on half densities, where P is, in the first place, a 
short range perturbation, that is, a first order self-adjoint differential operator with 
coefficients vanishing to second order at infinity. (Later, we show that there is a 
simple extension to metrics and perturbations of 'long range gravitational type', 
which includes the Newtonian or Coulomb potential and metrics of interest in 
general relativity.) We remark that the Riemannian half-density jcZgl^/^ trivializes 
the half-density bundle, and operators on functions can be regarded as operators 
on half-densities via this trivialization. Given A > 0, we solve for a kernel -R(A) on 
X"^ which satisfies 

(1.1) {H - X')R{\) = Km, 

where Km is the kernel of the identity operator on half densities. More precisely, 
we consider this equation on X^, which is the space X^ with the corner blown 
up. This allows us to use the scattering wavefront set at the 'front face' (the face 
resulting from blowing up the corner) to analyze singularities, which is an absolutely 
crucial part of the strategy. The kernel -R(A) is also required to satisfy a wavefront 
set condition at the front face, which is the analogue of the outgoing Sommerfeld 
radiation condition. 

We cannot find R{X) exactly in one step, so first we look for an approximation 
G{X) of it. The general strategy is to find G{X) which solves away the singularities 



of the right hand side. Km, of (1.1). Singularities should be understood both in 
the sense of interior singularities and oscillations, or growth, at the boundary, as 
measured by the scattering wavefront set. 

The first step is to find a pseudodifferential approximation which solves away 
the interior singularities of Km, which is a conormal distribution supported on the 
diagonal. This can be done and removes singularities except at the boundary of the 
diagonal, where H—X^ is not elliptic (in the sense of the boundary wavefront set). In 
fact, the singularities which remain lie on a Legendrian submanifold N* diag^, at the 



boundary of diagj, (see (4.5)). Singularities of G{X) can be expected to propagate 
in a Legendre submanifold i+(A) which is the bicharacteristic flowout from the 
intersection of N* diagj^ and the characteristic variety of iJ— A^. (The geometry here 
is precisely that of the fundamental solution of the wave operator in R"+^, which is 
captured by the intersecting Lagrangian calculus of Melrose-Uhlmann |0.) This 
Legendre has conic singularities at another Legendrian, L^(A), which is 'outgoing'. 
Thus, in view of the calculus of Legendre distributions of Mclrose-Zworski and the 
authors, the simplest one could hope for is that the resolvent on the real axis is the 
sum of a pseudodifferential term, an intersecting Legendre distribution associated 
to (iV* diagjj, _L+(A)) and a Legendre conic pair associated to (L+(A), U^{X)). This 
is the case: 



Theorem 1.1. Let H he a short range perturbation of a short range scattering 
metric on X . Then, for A > 0, the outgoing resolvent R{}? + iO) lies in the class 
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(4.4), that is, it is the sum of a scattering pseudodijjerential operator of order —2, 
an intersecting Legendre distribution of order —1/2 associated to {N* diagj,, L+(A)) 
and a Legendrian conic pair associated to (L-|_(A), L'^(A)) of orders —1/2 at L+(A), 
(n — 2)/2 at L'^{X) and {n — l)/2 at the left and right boundaries. 

If H is of long range gravitational type, then the same result holds except that the 
Legendre conic pair is multiplied by a complex power of the left and right boundary 
defining functions. 

This theorem was already proved in our previous work [Q , so it is the method that 
is of principal interest here. By comparison withM , the proof is conceptually much 
shorter; it does not use any results from ||l^ or [[l2| , though it makes substantial use 
of machinery from [p^ . But the main point we wish to emphasize is that the proof 
works directly on the spectrum and nowhere uses the limiting absorption principle, 
a method of attack that we think will be useful elsewhere in scattering theory. It 
seems that things which are easy to prove with this method are difficult with the 
limiting absorption principle, and vice versa. For example, it is immediate from our 
results that if / is compactly supported in the interior of X, then u = R{X^ + iQ)f 
is such that x~'^^^~^^^^e~^^/^u g C°°{X), while it is not so easy to see that the 
resolvent is a bounded operator from x^L^ to x^^L^ for any I > 1/2. Using the 
limiting absorption principle, it is the second statement that is much easier to 
derive (following [|l| for example). Thus, we hope that this type of approach will 
complement other standard methods in scattering theory. 

In the next section, we describe the machinery required, including the scattering 
calculus on manifolds with boundary, the scattering-fibred calculus on manifolds 
with codimension two corners, and Legendre distributions in these contexts. The 
b-double space, which is a blown up version of the double space Xy. X which carries 
the resolvent kernel, is also described. The discussion here is rather concise, but 
there are more leisurely treatments in |lj] and 

The third section gives a symbol calculus for Legendre distributions on manifolds 
with codimension two corners. This is a straightforward generalization from the 
codimension one case. 

The fourth section is the heart of the paper, where we construct the parametrix 
G(A) for the resolvent kernel. Propagation of singularity theorems show that the 
simplest space of functions in which one could hope to find the resolvent kernel is 



given by (4.4). We can in fact construct a parametrix for the resolvent in this class. 
In the fifth section this is extended to long range metrics and perturbations. 

In the final section we show that one can modify the parametrix so that the error 
term E{\) is such that Id-|-£'(A) is invertible. This is done by showing that the 
range of H—X^ on the sum of C°°(X) and G{X)C°"{X) is dense on suitable weighted 



Sobolev spaces. Thus the parametrix may be corrected to an exact solution of (1.1). 
Such a result also shows the absence of positive eigenvalues for H. Finally, we show 
that the kernel so constructed has an analytic continuation to the upper half plane 
and agrees with the resolvent there. 

Notation and conventions. On a compact manifold with boundary, X, we use 
C°°{X) to denote the class of smooth functions, all of whose derivatives vanish at 
the boundary, with the usual Frechet topology, and C~°°{X) to denote its topo- 
logical dual. On the radial compactification of M" these correspond to the space 
of Schwartz functions and tempered distributions, respectively. The Laplacian A 
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is taken to be positive. The space L'^{X) is taken with respect to the Riemannian 
density induced by the scattering metric g. This density has the form adxdy/x"^^ 
near the boundary, where a is smooth. 

Acknowledgements. We wish to thank Richard Mehose and Rafe Mazzeo for 
suggesting the problem, and for many very helpful conversations. A. H. is grateful 
to the Australian Research Council for financial support. A. V. thanks the NSF for 
partial support, NSF grant #DMS-99-70607. 



2. Preliminaries 

2.1. Scattering calculus. Let X be a manifold with boundary SAT = y. Near the 
boundary we will write local coordinates in the form (x, y) where a; is a boundary 
defining function and y are coordinates on Y extended to a collar neighbourhood 
of dX. 

We begin by giving the definition of a scattering metric. The precise requirements 
for the metric (and many other things besides) are easiest to formulate in terms 
of a compactification of the space. Taking the function x = as a boundary 
defining function and adding a copy of y at a; = yields a compact manifold. A, 
with boundary dX = Y. Then the definition of scattering metric is given in terms 



of X in Definition 2.1 below. Regularity statements for the metric coefficients are 
in terms of the C°° structure on A; this is a strong requirement, being equivalent to 
the existence of a complete asymptotic expansion, together with all derivatives, in 
inverse powers of r as r ^ oo. The benefit of such a strong requirement is that we 
get complete asymptotic expansions for the resolvent kernel, and correspondingly, 
mapping properties of the resolvent on spaces of functions with complete asymptotic 
expansions. 

Definition 2.1. A (short range) scattering metric on A' is a Riemannian metric g 
in the interior of X which takes the form 

, dx'^ h' 

(2.1) g= + 

x^ x^ 

where h' is a smooth symmetric 2-cotensor on A which restricts to the boundary to 
be a metric honY | [To[ |. A long range scattering metric is a metric in the interior 
of A which takes the form 

dx^ h' 

(2.2) 5 = aoo— + — , 

x^ X 

where aoo is smooth on A, aoo = 1 + 0{x), and h' is as above |p^ . If ooo = 
1 — cx + 0{x^) for some constant c we call g a gravitational long range scattering 
metric. 

Examples. Flat Euclidean space has a metric which in polar coordinates takes 
the form 

rfr^ + r^dui'^ , 

where duj'^ is the standard metric on 5*"^^. Compactifying Euclidean space as 
above, we obtain a ball with x = r^^ as boundary defining function, and then the 
flat metric becomes 

dx^ duj"^ 
x'^ ' 

which is a short range scattering metric. 



THE RESOLVENT FOR LAPLACE-TYPE OPERATORS 



5 



The Schwartzschild metric on R" takes the form near infinity 

1 - rfr^ +r^dtJ^ 

which under the same transformation leads to a gravitational long range scattering 
metric 

l-2mx — + — . 

X'^ 

The constant m = c/2 is interpreted as the mass in general relativity. 

The natural Lie Algebra corresponding to the class of scattering metrics on X 
is the scattering Lie Algebra 

Vsc(X) ^ {V \ V = xW, where W is&C°° vector field on X tangent to dX}. 

Clearly this Lie Algebra can be localized to any open set. In the interior of X, 
it consists of all smooth vector fields, while near the boundary it is equal to the 
C°°(A')-span of the vector fields x^dx and xdy^. Hence it is the space of smooth 
sections of a vector bundle, denoted ^^TX, the scattering tangent bundle. Any 
scattering metric turns out to be a smooth fibre metric on ^'^TX. The dual bun- 
dle, denoted ^'^T*X, is called the scattering cotangent bundle; near the boundary, 
smooth sections are generated over C°°{X) by dx/x'^ and dyi/x. A general point in 
^"^TpX can be thought of as the value of a differential d{f /x) at p, where / S C°°{X), 
and in terms of local coordinates {x, y) near dX can be written rdxjx^ + [iidyijx^ 
yielding local coordinates {x,y,T,iJ,) on ^'^T*X near dX. 

The scattering differential operators of order fc, denoted Diffg^(Ar), are those 
given by sums of products of at most k scattering vector fields. There are two 
symbol maps defined for P e Dif['^^{X). The first is the 'usual' symbol map, denoted 
here crf„t(P), which maps to S''=C'=r*X)/S''=-i(''=T*A:), where S^{'"=T*X) denotes 
the classical symbols of order k on ^'^T*X. The second is the boundary symbol, 
cra(P) G S''=("'=T^;^X), which is the full symbol of P restricted to x = 0. This is well 
defined since the Lie Algebra Vsc(-'^) has the property [Vsc(Ar), Vsc(A^)] C a;Vsc(-'^), 
so commutators of scattering vector fields vanish to an additional order at the 
boundary. Dividing the interior symbol a^^^{P) by \^\g, where | • |g is the metric on 
scrp*j^ determined by the scattering metric, we get a function on the sphere bundle 
of '''^T*X. This may be combined with the boundary symbol into a joint symbol, 
jsc{P)j ^ function on Csc{X) which is the topological space obtained by gluing 
together the sphere bundle of ^'^T*X with the the fibrewise radial compactification 
of ^'^Tg-^X along their common boundary. 

The scattering pseudodifferential operators are defined in terms of the behaviour 
of their Schwartz kernels on the scattering double space X^^, a blown up version of 
the double space X^. This is defined by 

Xl = [X^;{dXf] and 

Xi = [X^:ddmg^], 

and diagij is the lift of the diagonal to X^. The lift of diagj^ to X^^ is denoted 
diagg;,. The blowup notation [;] is that of Melrose: see (sj or |Q. The boundary 
hypersurfaces are labelled lb, rb, bf and sf; see figure |^. The scattering pseudo- 
differential operators of order k, acting on half densities, "^^^{X), are those given 

by KDs^c-valued distribution kernels which are classical conormal at diagg^ , of order 
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Figure 1 . The b-double space and the scattering scattering space 
k, uniformly to the boundary, and rapidly vanishing at lb, rb, bf. (Here KDgc is 



the puUback of the bundle 7rf"'=f]^(X) ® 7r;"<=fi^(X) over ^ X^.) The space 
is defined to be cc' "^^^{X). 
The joint symbol map j^^ extends from DiS^^{X) to ^^^{X) multiplicatively, 



q G Csc{X) if j^^{A){q) ^ 0. It is said to have elliptic interior symbol (boundary 
symbol) if j^^{A) does not vanish at fibre- infinity (spatial infinity), and is said to 
be totally elliptic if j^^{A) vanishes nowhere. The characteristic variety of A, S(A), 
is the zero set of i^c(^)- 

The scattering wavefront set of a tempered distribution u E C^°°{X) (the dual 
space of C°°(X), the space of smooth functions on X vanishing with all derivatives 
at the boundary) is the closed subset of Csc{X) whose complement is 



■'WF{uf = {qe Csc(X) I 3A e elliptic at q such that Au G C°°{X)}. 



The interior part of the wavefront set (at fibre-infinity) is a familiar object, the 
standard wavefront set introduced by Hormander (except that each ray of the stan- 
dard wavefront set is thought of here as a point in the cosphere bundle). In this 
paper we arc mostly interested in the part of the scattering wavefront set at spatial 
infinity. In fact, the operators H we shall study will have elliptic interior symbol, 
uniformly to the boundary, so in view of the next theorem, solutions of {H—a)u = 
must have wavefront set contained in the part of Csc{X) at spatial infinity, which 
we denote K (that is, K = "'^T^j^X). 

There is a natural contact structure on K induced by the symplectic form w on 
T*X. Writing lu in terms of the rcscaled cotangent variables r, and contracting 
with the vector x^dx yields the 1-form 

X = I'x'^d^,^ = dr + ^i- dy, 

which is nondegenerate, and therefore a contact form. A change of boundary defin- 
ing function x' = ax changes x by a factor a~^, so the contact structure is totally 



jt{A).jZ{B)=j^,+^{AB) 



and such that there is an exact sequence 




(2.4) 



THE RESOLVENT FOR LAPLACE-TYPE OPERATORS 7 

well-defined. Given a Hamiltonian h on K, the Hamiltonian vector field on K 
determined by the contact form is 

dh d / dh dh\ d dh \ d 

This is the same as x~^Vh restricted to a; = 0, where Vh is the Hamilton vector 
field on ^'^T*X induced by h. Integral curves of this vector field are called bichar- 
acteristics of h (or of A, if h is the boundary symbol of A). 

Under a coordinate change, the variables r and fj, change according to 

, , dy da , 

T — ar, ji — ail-—— — t— — - x = ax. 
ay ay 

Since a > 0, this shows that the subset 

(2.5) K^ = {{y,T,ij)^K\T<Q} 

is invariantly defined. This is important in the definition of the outgoing resolvent, 
see (13). 

The boundary part of the scattering wavefront set behaves very much as the 
interior wavefront set part behaves, and in particular we have a propagation of 
singularities result for operators of real principal type: 

Theorem 2.2 (Melrose). Suppose A e ^^^^{X) has elliptic interior symbol, and 
real boundary symbol. Then for u G C°°{X°) n C^°°{X), we have 

(2.6) WF,,(^u) C WF,e(w), 

WF,,(u)\WF,e(Au) CS(A), 

WFsc(m) \ WFsc(^m) is a union of 
maximally extended bicharacteristics of A inside S(A) \ WFsc(^u). 

Thus, if Au = 0, then '"^WF(u) C K and consists of a union of maximally 
extended bicharacteristics of A inside S(A). 

As well as a boundary principal symbol defined on K, scattering pseudodiffer- 
ential operators also have a boundary subprincipal symbol. This is the 0{x) term 
of the full symbol at the boundary when the operator is written in Weyl form. It 
is important to keep in mind that it depends on a choice of product structure at 
the boundary; it does not enjoy quite the same invariance properties as does the 
standard (interior) subprincipal symbol. A practical formula to use for differential 
operators with symbol in left-reduced form, ie, such that 

a{a{x,y){x'^Dxy{xDy)°')^a{x,y)T'li°', D = -id, a = (ai, . . . , a„_i), 

is that for a{P) — p{y, r, fi) +xq{y, r, /i) + 0(x^), the boundary subprincipal symbol 
of P is given by 

Lemma 2.3. Let g be a short range scattering metric, let x be a boundary defining 
function with respect to which g — dx^ / x^ + h' / x"^ , and let H be a short range 
perturbation of the Laplacian with respect to g. Then in local coordinates {x,y), the 
sub-principal symbol of H vanishes at fi — 0. 



(2.7) 

and 

(2.8) 
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Proof. The operator H may be written 

H = {x^Dxf + {n-l)ix^Dx + x^Ag + a,j{x,y)x^Dy^Dy^+Q, Q £ x^DiS^^iX). 
Thus the left-reduced symbol as above is 

(T(-ff) = + hijiiiiij + X ((n — l)ir + OijUifij) + 0{x'^). 
Hence the sub-principal symbol is 

dhi 

(2.10) cTsuhiH) = i-Q^fJij + aij^iHj, 

which vanishes when /i = 0. □ 

We now define the gravitational condition for the perturbation P. 

Definition 2.4. A first order scattering differential operator P on X is said to be 
short range if it lies in x"^ T)\S\^{X) , and long range if it lies in xT)iS\^{X). Let g 
be a scattering metric and x a boundary defining function with respect to which g 



takes the form (2.1) or (2.2). P is said to be of long range gravitational type with 



respect to g if it has the form 

P = X \y ai{xdy-) -\- hx^dx + c 




near x — Q, where a,;, h and c are in C°°{X), and for some constants bo and cq, 
h = ho + 0{x) and c = co + 0{x). 

The point of the short range condition is that then the subprincipal symbol of 
both H = A + P vanishes at the radial sets /x = 0, r = ±A oi H — A^, whilst in the 
long range gravitational case, the subprincipal symbol is constant. In the general 
long range case, the subprincipal symbol is an arbitrary function on the radial set, 
which causes some inconvenience (but not insuperable difficulties) in constructing 
the parametrix for {H — — iO)^^. 

2.2. Legendre distributions. An important special case that occurs often is that 
'^'^WF{u) is a Legendre submanifold, or union of Legendre submanifolds, of K; 
moreover, in many cases, m is a Legendre distribution, which means that it has a 
WKB-type expansion, the product of a oscillatory and smooth term, as discussed 
below, which makes it particularly amenable to analysis. 

We let dimX = n, so that dimK = 2n — 1. A Legendre submanifold of -fiT is a 
submanifold G of dimension n — 1 such that x I" G = 0. Such submanifolds have 
several nice properties. One is that if a Hamiltonian, h, is constant on G then its 
Hamilton vector field is tangent to G. Another is that Legendre submanifolds may 
be generated in the following way: If _F is a submanifold of dimension n — 2, such 
that X vanishes on F, and if the Hamilton vector field of h is nowhere tangent to F, 
then the union of bicharacteristics of h passing through F is (locally) a Legendre 
submanifold. 

Let G be a Legendre submanifold, and let q € G. A local (nondegenerate) 
parametrization of G near q is a function 4>{y, v) defined in a neighbourhood of 
ya & y and uq G M*"', such that dv<t> — Q &i q' = (yo,vo), q — {y,d(^x,y){4'/^)) at q', 
(j) satisfies the nondegeneracy hypothesis 

are linearly independent at G = {{y, v) \ dy(f> = 0}, 1 < i < fc. 
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and near q, 

(2.11) G = {{y,d^,,y)[^)\{y,v)eC}. 

A Legendre distribution of order m associated to G is a half-density of the form 
u = {uq + XljLi ^j)^) where v is b, smooth section of the scattering half density 
bundle, uq € C°°{X), and Uj is supported in a coordinate patch {x,y) near the 
boundary, with an expression 



m+n/i-k/2 



[ e''t'^^y'''^'''aj{x,y,v)dv, 



where locally parametrizes G and aj G C°°{X x K*^), with compact support in 

V. Melrose and Zworski showed that Uj can be written with respect to any local 
parametrization, up to an error in C°°{X). The set of such half-densities is denoted 
I"'iX,G;'"'ni). The scattering wavefront set of u e /'"(X, G; ^^fi^) is contained 
in G. 

An intersecting Legendre distribution is associated to a pair of Legendre sub- 
manifolds, L = (Lo,ii), where Li is a manifold with boundary such that Lq and 
Li intersect cleanly at dLi. A local parametrization of {Lq, Li) near q Cz LqH Li is 
a function (j){y, v, s) defined in a neighbourhood of q' = (t/o, uo, 0) in y x M'' x [0, oo) 
such that dy(f) = a,t q', q = {y,d(^x^y-^{(p/x)){q'), (p satisfies the nondegeneracy 
hypothesis 

/ d<j)\ . 
ds, dd), and d{ — — are linearly independent at g' , 1 < i < k, 

\aviJ 

and near 

Lo = {(y, )\s = 0,dy(l) = 0}, 

Li = { {y, d(^^y) (^)\s>Q,ds4> = Q, dy(j> = 0]. 

A Legendre distribution of order m associated to L is a half-density of the form u = 
Mo+(X]^i "j)^' where visa, smooth scattering half-density, uq G r^{X, Li\ ^'^VL^) + 

Lo; ^'^VL^) (the subscript c indicates that the microlocal support does not 
meet the boundary of ii), and Uj is supported in a coordinate patch {x,y) near 
the boundary, with an expression 



1)/^^ ds J e''t'i^y'''''^'''aj{x,y,v,s)dv, 



where (jij locally parametrizes {Lq.Li) and aj e C°^{X x R'^ x [0,00)), with com- 
pact support in v and s. Again, Uj can be written with respect to any local 
parametrization, up to an error in C°°(X). The set of such half-densities is denoted 
I"^{X, L; '^'^ria ). The scattering wavefront set of u S I"^{X, L; ^"^n^) is contained in 

A Legendre distribution associated to a conic Legendrian pair is associated to a 
pair of Legendre submanifolds G = {G,G^) where G" is a projectable Legendrian 
(that is, the projection from ^'^T*X to F is a diffeomorphism restricted to G") and 
G is an open Legendrian submanifold such that G \ G is contained in G** and G 
has at most a conic singularity at G'. We further assume that r ^ on G", so 
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that we may change coordinates to a new boundary defining function such that 
is parametrized by the phase function 1. In these coordinates, the condition that 
G has a conic singularity at G" means that G hfts to a smooth submanifold with 
boundary, G, on the blown-up space 

(2.12) r"; {^ = 0,M-0}], 



intersecting the front face of ( p.l2| ) transversally. In local coordinates (a;, y, r, /i), 
coordinates near the front face are 

(2.13) x/l/il, T, and /i, 

and we require that G is given by the vanishing of n smooth functions of these 
variables with linearly independent differentials, and that dj/ij 7^ at dG. 

A local parametrization of G near g <E G n G* is a function (j){y,v,s) = 1 + 
sijj{y,v,s) defined in a neighbourhood of q' = {yo,vo,0) in F x M'^ x [0,oo) such 
that (j>o parametrizes G" near q, dycj) — &i q\ q — {y,d(^j. y^((j)/x)){q'), (j> satisfies 
the nondegeneracy hypothesis 

/ dip \ 

ds, dip, and d\ — — ) are linearly independent at g', 1 < i < k, 
\dviJ 

and near g, 

G = {(0, y, -0, sdyTp, dyip) \ dyCp = 0, dsip = 0}, 
in the coordinates ( 2.13| ). A Legendre distribution of order {m,p) associated to 



(G, G") is a half-density of the form u = uq + (X^iLi^i)^; where 1/ is as above, 
uo e IJJ'iX, G; '"'ni) + IP{X, G»; "'^f^s) (the subscript c indicates that the microfocal 
support does not meet G**), and uj is supported in a coordinate patch {x,y) near 
the boundary, with an expression 

°° r , ^, / r \ 'n+n/4-(fe+l)/2 

d5 / e*0.(y,-.^)/-a,(2/,i;,a;/s,s)(^) sP+^^Z^-^ dv. 







where (pj locally parametrizes (G, G*) and aj S C°°{X x R'^ x [0, 00) x [0, 00)), with 
compact support in x/ s and s. Here Uj can be written with respect to any local 
parametrization, up to an error in IP{X,G'^]^^Vt^). The set of such half-densities 
is denoted I'^-p{X, G; '^'^f^^). The wavefront set of u e I'^iX, G; """"n^) is contained 
in GU G«. 



2.3. Codimension 2 corners. In this subsection we briefly review the extension 
of the theory of Legendre distributions to manifolds with codimension 2 corners 
and fibred boundaries given in jj]. 

Let M be a compact manifold with codimension 2 corners. The boundary hyper- 
surfaces will be labelled mf , Hi, ... , H4, where the Hi are endowed with fibrations 
TTj : Hi — > Zi to certain closed manifolds Zi and mf (the 'main face') is given the 
trivial fibration id : mf mf. The collection of fibrations is denoted It is 
assumed that HiO Hj = % \i i ^ j. It is also assumed that the fibres of tt^ intersect 
Hi n mf transversally and therefore induce a fibration from Hi n mf — > Zi. Further, 
it is assumed that a total boundary defining function x is given, which is distin- 
guished up to multiplication by positive functions which are constant on the fibres 
of dM. 
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Near H D mf , where H = Hi for some i, there are coordinates a;i, X2, 2/1,2/2 such 
that xi is a boundary defining function for H, X2 is a boundary defining function 
for mf, X1X2 = X, and the fibration on H takes the form 

(yi, 2:2,2/2) 1-^ 2/1- 

Associated with this structure is a Lie Algebra of vector fields 

Vs4.(M) = {V \ V gC°^,V is tangent to $ at dM, V{x) = 0{x^)}. 

This is the space of smooth sections of a vector bundle, denoted ^"^TM. The dual 
space is denoted ^^T*M. A point in ^^T*M may be thought of as a differential 

d{f /x) at p, where / is a smooth function on M constant on the fibres at dM. A 
basis for ^"^TpM, for p e M near mfflH, is given by dx/x"^, dxx/x, dyi/x, dy2/x^. 
Writing q G ^■^T*M as 

dx dxi dyi dy2 

q = T—^ + n h fll 1- /X2 

X-' X X X2 

gives coordinates 

(2.14) (xi, X2, 2/1, 2/2, r, n, /ii, /X2) 

on '^*r*M near mf nif. 

The differential operators of order at most k generated over C°°(M) by Vs$(M) 
are denoted Diffg^(M). Near the interior of mf, the Lie Algebra Vs$(M) localizes 
to the scattering Lie Algebra Vsc{M), where M denotes the noncompact manifold 
with boundary M \ (JiHi. Consequently, we have a boundary symbol (ts[P), P € 
Diff^$(M) taking values in S''{'"'T*^M) over the interior of mf. In fact the symbol 
extends to an element of S''{^^T^fM) continuous up to the boundary of mf. 

For each fibre F of H, there is a subbiuidlc of ^'^TM consisting of all vector 
fields vanishing at F. The annihilator subbundle of "*T*M is denoted "'=T*(if; F) 
since it is isomorphic to the cotangent space of the fibre. The quotient bundle, 
^^T*M/'"'T*{H\ F) is denoted "*7V*Zj since it is the pullback of a bundle over Zj. 
The fibration tTj induces a fibration 

(2/1, 2/2, T,Ti, ^1,^2) ^ (r,yi,/zi). 

We next describe three types of contact structures associated with the structure 
of M. Since Vs$(M) is locally the scattering structure near the interior of mf, 
there is an induced contact structure on ^*T*M over the interior of mf. In local 
coordinates, the contact form looks like 

X = '-3,29^ [lo) = dr + Tidxi + ^idyi + xi/i2C?2/2. 

We see from this that at cci = 0, x is degenerate. However, restricted to mf flff, x 
is the lift of a form xz, on ^^N*{Zi), namely dr + Hidyi, which is nondegenerate 
on '^'^N*{Zi). This determines our second type of contact structure (one for each 
i). The third type of contact structure is that on ^'^Tgp{Hi; F) induced by Vs<i>(M) 
for each fibre F of Hi, since it restricts to the scattering vector fields on each fibre. 
In local coodinates, this looks like dri + /i2C?2/2. 

Using these three contact structures we define Legendre submanifolds and Le- 
gendre distributions. 
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Definition 2.5. A Legendre submanifold G of "*T*M is a Legendre submanifold 
of '^^Tj^fM which is transversal to ^^T^^nmf^^ '^^'^h which the map ( ^.15 ) 

induces a fibration from dG to Gi, where Gi is a Legendre submanifold of ^'^N*Zi, 
whose fibers are Legendre submanifolds of ^"^TgpF. 

A projectable Legendrian (one such that the projection from ^'^T^fM mi is 
a diffeomorphism when restricted to G) is always of the form 

graph (d(^))={(f,4^))|l7emf} 

for some smooth function (p constant on the fibres of dM. We then say that (f> 
parametrizes G. In general, let G be a Legendre submanifold of ^'^T*AI, and let 
q ^ G. li q lies above the interior of mf, then a local parametrization of G near q is 
as described in the previous subsection, so consider q G dG lying in ^^T^^M, where 
Ijq G mf nH. A local (nondegenerate) parametrization of G near g is a function 
4>{xi, j/i, 2/2, V, w) of the form 

(2.16) (I){xi,yi,y2,v,w) = (l)i{yi,v) + Xi(l)2{xi,yi,y2,v,w), 

defined in a neighbourhood of q' — {Ijq, vq^wo) G mf xM'''i+'^^ , such that dy^^cf) = 
at q' , 

q={y^,d{^))e.iq' 

in local coordinates ( ^.14 ) , satisfies the nondegeneracy hypothesis at q' 
(2.17) 

d{yi,v)^^ij = 1, ■ ■ ■ , and d(v2,w)^^^-, / = !,■■■ linearly independent, 
and near g, 

(2.18) G = {(y,d(^)) |40 = d^^ = O}. 

A Legendre distribution of order (m; ri, . . . r^) associated to G is a half-density such 
that for any Vi E C°°(M) whose support does not intersect Hk, for k ^ i, ViU is 
of the form u = Uq + (X^jLi + Sjli ^j)^' where v \s a smooth section of the 
half-density bundle ""^n^ induced by "*r*M, uq € C°°(X; ''*r25), and u,, u'^ have 
expressions 



u,(xi,a;2,yi,2;2) = / e''^^("i*'«^'"''")/"a,(a;i, a;2, yi, 1/2, i;, H 

(2.19) ^ 

^™-(fe+fe')/2+A^/4^n-fc/2+A./4-/./2^^^^ 

with N = dimM, G C^([0, e) x [/ x R'''+'"''), f7 open in mf, the dimension of 
the fibres of Hi and (pj a phase function parametrizing a Legendrian G on U , and 

(2.20) u'(a;i,yi,z) = / e"^^(«i'"')/^aj(a;, yi, z, u;)a;'^-~'=/2+^/''-/-/2 du; 



with TV = dimM, Oj G C^([0,e) x U x R*-'), C/ open in i7, fi as above, a phase 
function parametrizing the Legendrian Gi. 
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Definition 2.6. A Legendre pair with conic points, (G, G"), in ^'^T*M consists 
of two Legendre submanifolds G and G' of ^'^T*M which form an intersecting 
pair with conic points in ^'^T*-^M such that for each Hi the fibrations of G and 
G" induced by (2.15) have the same Legendre submanifold Gi of ^'^N* Zi as base 
and for which the fibres are intersecting pairs of Legendre submanifolds with conic 
points of "'=r^;^F. 

The Legendrian G« is required to be projectable, so it parametrized by a phase 
function 4>{jj) which is constant on the fibres of dM. Thus, x' — x/cf) is another 
admissible total boundary defining function. With respect to x' , G" is parametrized 
by the function 1. Thus, without loss of generality we may assume that coordinates 
have been chosen so that G' is parametrized by 1. This simplifies the coordinate 
form of the blowup (2.12). Coordinates near dG then are 

(2.21) y = {xi,yi,y2), t, Ti/|/i2|, fJ.i/\fJ.2\, h ^ fJ-2/\m\, and |^2|, 

the last of which is a boundary defining function for G (see [||). 

As a consequence of definition 2.6, G is a compact manifold with corners in 

(2.22) {x = 0, = /i2 = n= 0}], 

with one boundary hypersurface at the intersection of G and ^'^T^^^jj AI for each 
i (which are mutually nonintersecting, since the Hi are mutually nonintersecting) , 
and one at the intersection of G and the front face of the blowup in ( 2.22 ). If q 
lies in the interior of G then the situation is as for Legendrians in the scattering 
setup. If q is on the boundary of G, but does not lie over Hi for some i, then the 
situation is as for Legendrian conic pairs as in the previous subsection. If q is on 
the boundary of G but not in G'* then the situation is as above. Thus the only 
situation left to describe is if q is in the corner of G, lying above ^ ™f ^^^i say. 

A local parametrization of (G, G**) near q (in coordinates as chosen above) is a 
function 

(I>{^i,yi,y2,s,'w) = 1 + sxiip{xi,yi,y2,s,w), 
with defined in a neighbourhood of q' = {y^, 0, wq) C M x [0, oo) x M''', such that 
dw'ip = at q' , 

f dip\ , 
dy^t/j and (i(j,2,^) I ^ — -I are linearly independent at q , 

and such that near q G, 

A rr- , dxAxiip) dy^{xiip) -r—- 

\dy^'lp\ \dy^i^\ 

in the coordinates ( 2.21 ). 

A Legendre distribution of order {m,p; ri, . . . r^) associated to (G, G") is a half- 
density such that for any Vi G C°°{M) whose support does not intersect Hk, for 
k i, ViU is of the form m = uq + (X)j=i + "j")^' where uq G C°°(X; '^^ils)^ 

and Uj, u'j have expressions 



(2.23) 



Uj{xi,X2,yi,y2)^ ds j du; e^"*^(^^'^^'^^''^'"')/"a,(xi, s, xs/s, z/i, ?/2, ^i^) 

,^.m-(k' + l)/2+N/4 ^,_fe/2+Ar/4-/./2 
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with = diniM, G C^([0, e) x [/ x U open in mf, fi the dimension of 

the fibres of Hi and (f> j a phase function parametrizing a Legendrian G on [/, and 



where w' is as in (2.2C) 



2.4. The b-double space. Here we analyze the b-double space X^, where X is a 
compact manifold with boundary, from the perspective of manifolds with corners 
with fibred boundaries. The manifold with corners has three boundary hyper- 
surfaces: lb and rb, which are the lifts of the left and right boundaries dX x X, 
X X dX of X'^ to X^, and bf, coming from the blowup of (dX)^ (see figure 0). 
Thus, lb and rb have natural projections to dX. The fibres of lb and rb meet bf 
transversally, so we may identify bf as the 'main face' mf of X^. Given coordinates 
{x,y) or z on X, we denote the hft to X^ via the left, resp. right projection by 
{x',y') or z', resp. {x",y") or z" . We may take the distinguished total boundary 
defining function to be x' , for a = x'/x" < C and x" , for a > C~^. These are 
compatible since their ratio is constant on fibres on the overlap region < a < C 
(this is trivially true since the fibres of bf are points). 

These data give X^ the structure of a manifold with corners with fibred boundary 
as defined above. The s$- vector fields then are the same as the sum of the scattering 
Lie Algebra Vb(X) lifted to X^ from the left and right factors. 

On X^, and ^'^T*X^, it is most convenient to use coordinates lifted from X and 
^'^T*X. Near lb, but away from bf, we use coordinates {x',y',z") and coordinates 
(r', /i', C") on «*T*X2 where we write a covector q e ^■^T*X^ as 

q^r'%+,'.'4 + C"-dz". 
x' x' 

Similarly near rb, but away from bf, we use coordinates (z', a;", y"; C', t", m")- Near 
lb n bf , we use (x", ct, y', y") with corresponding coordinates (r, k, /z', by writing 

q e "*r*x^ as 

dx' da I dy' „ dy" 

x' x' x' x" 

However, we may also use scattering cotangent coordinates (r', r", /x', /.t") lifted 
from ^'^T*X, where we write 

,dx' ,,dx" , dy' ,, dy" 
x'^ x"^ ^ x' ^ x" 

This gives 

(2.24) r' = T + CTK r" = -K. 

The coordinates {x" ,a,y' ,y") hold good near bf as long as we stay away from rb, 
when we need to switch to {x' , ,y' ,y"). The cotangent coordinates (r', r", //', /i") 
are good coordinates globally near bf; notice that the roles of {t,Ti, ^i, ^2) are 
played by {t' , t" , fi' , fi") near lb and (r", r', /x", near rb. 

The operator H can act on half-densities on X^ by acting either on the left or 
the right factor of X; these operators are denoted Hi and H^ respectively. For 
H = A + P, where P £ xDiSl^{X), the Hamilton vector field induced by Hi and 
the contact structure on ^'^T^fX^, with respect to x' , takes the form 

d d d f dh' d dh' d \ 

(2^25, V, ^ 2^4 + 2rV A - + (|t_| - „■ . 
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Similarly, the Hamilton vector field induced by and the contact structure on 
''^T^fX^, with respect to x", takes the form 

(2.26) 



Notice that Vi and Vr commute. 



dh" d dh" d 



dfi" dy" dy" dfi" 



h" = h{y",f,"). 



3. Symbol calculus for Legendre distributions 

3.1. Manifolds with boundary. Let X be a manifold with boundary of dimen- 
sion TV, and let 



,-fc/2-n/4 



(y'")/^a(a;,?/,w)dv 



dxdy 



G /™(X,G;"'=17: 



be a Legendre distribution of order to. Let C — {{y,v) \ dycj) — 0} and let A be 
a set of functions in (y,w)-space such that {\,dy(f) form local coordinates near C. 
We temporarily define the symbol relative to the coordinate system Z = {x, y) and 
the parametrization to be the half density on G given by 



(3.1) 



a^,^{u)^{a[Q,y,v)\C) 



d(dy(t>, A) 



d{y,v) 



Here we have used the correspondence ( |2.1l| ) between C and G. 
If we change coordinate system, the symbol changes by 



X 
X 



da 



(3.2) 
where 

If the parametrization is changed, then by Q, the symbol changes by 
(3.4) = CTS',^(u)e*'^('''s"''-^-'^'sn<,„,0)/4. 

the exponential is a locally constant function. We use these transformation fac- 
tors to define two line bundles, the i?-bundle over ^'^T*X which is defined by the 
transition functions (3.3), and the Maslov bundle over G which is defined by the 
transition functions ( [3.4| ) . (These bundles will be described in much more detail in 
§.) Defining the bundle S'["1(G') = \N*dX\"^-'^/'^ E (g) M{G) over G, we obtain 
an invariant symbol map from (|3.l|) 



cr™ : /™(X,G;'^'=17^) ->C°°(G;175 ® S'^'^HG)). 

The elements of the symbol calculus for Legendre distributions on manifolds with 
boundary have been given by Melrose and Zworski p^ : 

Proposition 3.1. The symbol map induces an exact sequence 

^ /"+i(X, G; """"n^) /™(X, G; """"n^) -> C°°(G, ® S'['"1(G)) -> 0. 

//p e 4»fe(X;'''=r23) andw e /'"(X,G;'''=175), then Pu e I'"'{X,G;'"'ni) and 

cr™(Pu) = (cr(P) \ G)f7™(w). 
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Thus, if the symbol of P vanishes on G, then Pu G /'"^^(X, G; '^'^O^). The symbol 
of order m + I of Pu in this case is 

t 1 N dp \ 

(3.5) - iCH, " + ~ 1") a7 + Psuh)cr"'{u) (E) \dx\, 

where Hp is the Hamilton vector field of p, the principal symbol of P, and Psuh is 
the subprincipal symbol of P. 

The symbol calculus for intersecting Legendre distributions is easily deduced 
from Melrose and Uhlmann's calculus of intersecting Lagrangian distributions. The 
symbol takes values in a bundle over LqULi. Let pi be a boundary defining function 
for dLi as a submanifold of Lo, and po be a boundary defining function for dLi as 
a submanifold of Li. To define the symbol, note that the symbol on Lq is defined 
by continuity from distributions microsupported away from Li, and takes values in 

(3.6) 

pr'C°°(r!i/2(io) ® 5["+i/2l(Lo)) = Pi'^^C^inl^^Lo \ dLi) ® 5[™+i/2l(Lo)), 

while the symbol on ii defined by continuity from distributions microsupported 
away from dLi takes values in 

C°°(r!i/2(Li) ® ^ pl/^C°^{nl^\Li) ® 5['"l(Li)). 

Melrose and Uhlmann showed that the Maslov factors were canonically isomorphic 
on I/gnLi, so 5[™l(Lo) is naturally isomorphic to S^'"^\Li) over LoHii. Canonical 
restriction of the half-density factors to Lq D Li gives terms in C°° {il^ (Lq n Li) (g) 
S'^'^\Li)(g\Nl^dLi\-^/'^®\N*dX\^/'^ eindC'^{^l^LonLi)(gS^'^^{Li)(g\Nl^dLi\^/^ 
respectively. In fact \Nl^dLi\ (g) \Nl^dLi\ (g) \N*dX\~'^ is canonically trivial; an 
explicit trivialization is given by 

(3.7) {dpo,dpi,x~'^) 1-^ x^'^uj{Vp„,Vp,) ILqDLi, 

where V^. are the Hamilton vector fields of the functions pi extended into ^'^T*X, 
and Lu is the standard symplectic form. Thus the two bundles are naturally isomor- 
phic over the intersection. We define the bundle S^'^^L) to be that bundle such 
that smooth sections of Vll^"^ (L) iSi S^"^\L) are precisely those pairs (a, b) of sections 
of p^^C°°{n^/^{Lo) (g) S'[™+i/2](Lo)) and pl^^C°° {^1^^ {Li) (g ^^(L^)) such that 

(3.8) p^h = e*"/4(2^)i/ Vo'/'a at Lq D Li 

under the identification (3.7). The symbol maps of order m on Li and m + 1/2 on 
Lq then extend in a natural way to a symbol map of order m on L taking values in 
nl^\L)g,S^"'\L). 

Proposition 3.2. The symbol map on L yields an exact sequence 

(3.9) 0-^ /"+i(X,L;^'=f75) ^ L"\X,L;'"'ni) C°°{L,nl g, S^""^) 0. 
Moreover, if we consider just the symbol map to Li, there is an exact sequence 

(3.10) ^ r'+'^{X,L;'"'ni) + I"'+i{X,LQ;'"'ni) I"'{X,L;^''ni) 

^C°°(ii,f]5(g,S'M)^0. 
//p e 5rfe(X;'''=r23) andue I"'{X,L;'"'ni), then Pu e P''{X, L^'^'ni) and 

a™(Fw) = (cr(P) \ Ly"\u). 
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Thus, if the symbol of P vanishes on Li, then Pu is an element of P"^^ {X , L; ^^^^2 ) + 
P^{X,Lq;^^Q,^). The symbol of order m + 1 of Pu on Li in this case is given by 

For a conic pair of Legendre submanifolds G — (G, G"), with G the desingularized 
submanifold obtained by blowing up G', the symbol is defined by continuity from 
the regular part of G. The symbol calculus then takes the form 

Proposition 3.3. Let s he a boundary defining function for G. Then there is an 
exact seguence 

^ I"'+^-P[X,G\'"'n'^) G; ""Tl'^) s™-PC°°(G,^2j (»S'["1(G)) ^ 0. 

If P (^^^{X^'^'^Vt^), andue I"^P{X,G;'"'n^), then Pu e I"^P{X,G;'"'n^) and 

a"'{Pu) = (a{P) \ G)a"(w). 

// the symbol of P vanishes on G, then Pu G G; ^''^fis). The symbol of 

order to + 1 of Pu in this case is given by (|3.5| ). 

3.2. Codimension two corners. When we have codimension two corners, then 
essentially the same results hold by continuity from the main face. The symbol 
is defined as a half-density on G by continuity from the interior of mf, where the 
scattering situation applies. We must restrict to differential operators, however, 
since pseudodifferential operators have not been defined in this context. 

Let M be a manifold with codimension 2 corners with fibred boundaries, let 
N = dimM, and let G be a Legendre distribution. Let pi be a boundary defining 
function for Hi. The Maslov bundle M and the E-bundle are defined via the 
scattering structure over the interior of G and extend to smooth bundles over the 
whole of G (that is, they are smooth up to each boundary of G at ^"^T^.^^fM). 
Let 5[™1(G) = M{G)®E(g) |A^* mf l™"^/"* ^ \N* Hi\^~^ ® ■ ■ ■ ® \N* Hd^^-^ ''^ . 
Finally let r stand for (ri, . . . , r^), and let = Hi • 

Proposition 3.4. There is an exact sequence 

0^ /™+^''"(M, G;"*f75) /'"''■(Af,G;''*f75) -> p'"-'"C~(G, ^ij ® S'["1(G)) ^ 0. 
//P e Diff(Af;"*rj5) andu e /'"'""(Af,G;"*17^), then Pu & I"^^"^ {M^Gf^Vl^) and 

a"'{Pu)= (a(P) fG)a"(w). 

Thus, if the symbol of P vanishes on G, then Pu G I"'+^'''{M,G;'"^ni). The 
symbol of order m + 1 of Pu in this case is given by (|3.5| ) . 

For a conic pair of Legendre submanifolds G — (G, G'*), with G the desingularized 
submanifold obtained by blowing up G'*, the symbol calculus takes the form 

Proposition 3.5. Let s be a boundary defining function for G at G O GK Then 
there is an exact sequence 

(3.11) 0-^ /'"+i'P'''(Af,G;''*f]5) ^ /'".P;'-(M,G;''*f)5) 

^ p™-r^™~PC~(G,r!j (g> 5["1(G)) ^ 0. 
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//P e Diff(A/;"*r2^), and u £ /"'^^""(M, G; then Pu G /'"■^'""(M, G; 

and 

//i/ie si/m6o; o/P vanishes on G, then Pu e I"'+^'P'''{M,G;^'^n^). The symbol of 



order m + I of Pu in this case is given by (3.5) 



The proofs of these propositions are omitted, since they are easily deduced from 
the codimension one case. 

4. Parametrix construction 

In this section, we consider self-adjoint operators H of the form A + P, where 
A is the positive Laplacian with respect to a short-range metric on a compact 
manifold with boundary, X, and P £ Diffg^(X) is a short-range perturbation of 
A. In the following section, we consider metrics and perturbations of long range 
gravitational type. Let R{a) denote the resolvent {H — a)^^ of H. 



In this section, we directly construct a parametrix G(A) for (1.1) whose error 
term E{X) — {H — A^)G(A) — Id is compact. Using Fredholm theory and a unique 
continuation theorem we solve away the error, giving us a Schwartz kernel R{X). 
We then show that R{y/a) has an analytic continuation (as a distribution on X'^) 
to the upper half a plane which agrees with the resolvent i?(cr) there. This proves 
that R{X) and R{X^ + iO) coincide on the real axis. 

The distribution R{X) has the defining property that 

(4.1) {H - A2)P(A) = Id, 
as an operator on C°°{X;^'^ili), and that 

(4.2) ""WF'^iRiX)) C if- as defined in (O). 



Here X^ ^ X^\ {lb U rb} is regarded as an open manifold with boundary, so that 
we can talk about the scattering wavefront set over the interior of bf. Equation 



(4.2) is the microlocal version of the outgoing Sommerfeld radiation condition. (For 



example, if A > 0, e*"**/^ has wavefront set in K^, while e ^'^^^ does not.) 



Equation (4.1) means that the kernel of P(A), which we also denote by -R(A) by 



an abuse of notation, satisfies 

(4.3) [Hi - X')R{X) ^ Km, 

where Kid is the kernel of the identity operator, i.e. it is a delta distribution on the 
diagonal, and Hi is the operator H acting on the left factor of X in X x X. 

There are four main steps in the construction. First we find an approximation 
to -R(A) in the scattering calculus, Gi(A) G \1/~^(X), which removes the singularity 



on the diagonal in (4.3). This leaves an error which, when viewed on the b-double 
space X^, is singular at the boundary of the diagonal 9diagi,. In fact, it is Leg- 
endrian at a Legendre submanifold lying over d diagj, which we denote N* diag]-, 
(see ([4.5|)). We solve this error away locally near Sdiagj, using an intersecting 
Legendrian construction which is due (in the Lagrangian setting) to Melrose and 
Uhlmann ||ll|; the singularities inside N* dia,gy^r]Y.{Hi — A^) propagate in a Le- 
gendre submanifold Z/+(A). This Legendre submanifold intersects both lb and rb, 
and an 'outgoing' Legendre submanifold L*(A); (L+(A), L'^(A)) form an conic pair 
of Legendre submanifolds and we can find a conic Legendre pair which solves away 
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the error up to an error term which is Legendrian only at L'^(A), ie we can solve 
away the errors at i+(A) completely. Finally, this outgoing error is solved away, 
using a very standard argument in scattering theory, at lb and bf, leaving an error 
E{X) which is compact on weighted spaces x''L'^{X) for all I > 1/2. 
Thus, wc seek G(A) (and R{X)) in the class 

(4.4) 

*2.0(X)+/-5(iV* diag^,L+{X);^'n^) + r-2'^'^''^{L+{X),L*{X);'^ni) 

where the second term is an intersecting Legendrian distribution and the third is a 
Legendre conic pair with orders -1/2 at L+(A), (n- 2)/2 at L*{X) and (n - l)/2 



at lb and rb. In this class of distributions there is a unique solution R{X) to (4.1) 



and (4.2) 



To avoid cumbersome notation, Q will denote a generic correction to the para- 
metrix constructed so far, and E will denote a generic error. The values of these 
symbols is allowed to change from line to line. 

4.1. Pseudodifferential approximation. The first step in constructing G(A) is 
a very standard argument. We seek G'i(A) e ^sc^(X) such that 

{H^X^)Gi{X)^ld+E,{X), Si (A) e %,°"{X). 

This will mean that the error term -E'i(A) has a smooth kernel (times the standard 
half-density) on X^^, so that we have solved away completely the singularity along 
the diagonal. 

The standard elliptic argument applies here since the interior symbol of H — X^ 
is elliptic. Thus, we first choose any Q G 'i'^^{X) whose interior symbol is | • |g ^ = 
(a'^iH - X^))-\ Then 

{H - X^)Q = ld+E, Ee%^\X). 

Multiplying Q by a finite Neumann series (Id +E + • • • + E^^^) thus gives an error 
E^ G ^'g^'^(X). Taking an asymptotic limit gives us a Gi(A) e %^^{X) with the 
desired error term. 

4.2. Intersecting Legendrian construction. In the next step of the construc- 
tion we move to X^, and view the error -E'i(A) from the first step of the construction 
on X^ rather than X^^. On X^ it has a smooth kernel except at Sdiagj^ where it 
has a conic singularity. That is, at Sdiag^, = {x' — 0,a = x'/x" = l,y' ~ y"}, 
the kernel is a smooth (and compactly supported) function of a:', S* = (ct — l)/x', 
Y = {y' — y")/x' and y'; this is easy to see since these are smooth coordinates on 
sf C Xg^. Using the Fourier transform, we write 

^;i(A) = (^j e'^^y'-y"^-''+^''-^^*'>/''' a{x\y' ,T],t) dr]d?j v 

where a is smooth in all variables, and in addition Schwartz in (?7,i). The phase 
function (y' — y") ■ rj + (a — l)t parametrizes the Legendrian 

(4.5) N* diagb = {y' = y", cr = 1, ^' = -fi", t' = -t"}. 

Therefore, £'i(A) is a Legendre distribution of order associated to N* diagj^. (To 
be pedantic, £'i(A) does not fall strictly in the class of Legendre distributions as 



defined by Melrose and Zworski since its microsupport is not compact; from (4.5) 



we see that the microsupport is a vector bundle over Sdiag,-,. It is instead an 
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'extended Legendre distribution' as defined in Q. However this is of no significance 
since the symbol is rapidly decreasing in each fibre of the vector bundle, hence all 
constructions we wish to perform here are valid in this context.) 

Observe that ao{Hi — A^) = t'" + h{y',fi') — vanishes on a codimension one 
submanifold of N* diagj,, and does so simply. Consider the vector field V/ which 
is given by (2.25). Since r'^ + ft, = A^ 7^ on — A^), at least one of the 

coefficients of da and dr' in (2.25) is nonzero, so Vi is transverse to A^* diagj, at 
the intersection with I](i?; — A^). We define L°{X) to be the flowout Legendrian 
from TV* diag^, nl](77; — A^) with respect to V/, and i±(A) to be the flowout in 
the positive, resp. negative direction with respect to V/. Thus, at least locally 
near N* diag^,, i±(A) are smooth manifolds with boundary. Notice that by (2.24), 
N* diagi, is contained in t = and Vi{t) < at N* diag^,. Thus, at least locally 
near iV* diagj,, L'^{X} is contained in = {t < 0}. The global properties of 
Lj-(A) are studied in the next section; in this section we only work microlocally 
near N* diagi,. 

The first step in solving away the error -Ei(A) from the previous step is to find 
an intersecting Legendrian Q e /^^/^(X^, {N* diag^, L+(A)); ''*ri2 ) such that 

(4.6) 

(Hi - \')Q ~ E,iX) e I' {Xl , N* diagb ■^^n^)+li{Xl (N* diag^ , L+ ( A) ) , ^* 17^ ) , 

microlocally near N* diagi,. To do this we choose Q with symbol on N* diagj, equal 
to a9){Hi - A2)-icrO(£:^^^(A)). This is an admissible symbol on N* diag^ by (Q, 
and (|3.9|), since (Tq{Hi — A^) is a boundary defining function for L^{X) on N* diag,^. 



It determines the value of the symbol on i+(A) at the boundary by ( |3.8| ). We 
extend this symbol by requiring that the transport equation, (3.5), be satisfied. 
This equation is a first order linear ODE with smooth coefficients, so there is a 
unique solution in a neighbourhood of N* diag^,. Then the symbol of order —1/2 
of {Hi — A^)i?o ~ £'^^^(A) vanishes, and there is an additional order of vanishing on 
L+(A) since the transport equation is satisfied. Thus by (3.10) the error term is as 
in(|3)- 

We now show inductively that we can solve away an error Ek which is in 
jf'iN* diagb)+/'=+i/2(A^* diagb, L+(A)) with a term Qk e l''+'^/^{N* diag^, i+(A)), 
up to an error which is in /'^+^(iV* diag^,) + 7'=+3/2(_/V* diag^, i+(A)). The argu- 
ment is the same as above: we take the symbol of order k on N* diagj, equal to 
(y'^iHi — \^)~^ (7^ {Ek) ^ and the symbol on i+(A) to solve away the symbol of order 
A; + 1/2 of Ek when the transport operator is applied to it. Taking an asymptotic 
sum of Q and the QkS gives us an error term which is microlocally trivial near 
N* diagj,. By cutting off away from Sdiag^,, we obtain an error 

i?2(A) e ll'\L+{\)), 

where the subscript c indicates that the microlocal support is compact and disjoint 
from the intersection with N* diagj^. 

4.3. Structure of L{\). In this section we analyze the global structure of L°{X). 
This was defined as the flowout from TV* diagjj C\Y.{Hi — A^) by the vector field V/. 
In fact, it is quite easy to see that A^* diagi, nl](i?; — A'^) = A^* diag,^ nE(i/r — 
A^). Moreover, neither Vi nor Vr is tangent to A^* diagj, at any point contained 
in ll{Hi — A^), but the difference V/ — K is tangent to A'* diagj,. Since V/ and K 
commute, this shows that the flowout with respect to Vi is the same as the flowout 



THE RESOLVENT FOR LAPLACE-TYPE OPERATORS 



21 



with respect to Vr- We will soon see that the symbols of our parametrix on _L+(A), 
defined so as to satisfy the left transport equation, also satisfy the right transport 
equation. 

It is convenient to write down L°{X) explicitly. Indeed, the computation of 
Melrose and Zworski can be applied with a minor change (that takes care of the 
behavior in a) to deduce that 

(4.7) 

L°iX)^{ie,y',y",T',T",^l',^i")■■ 3{y,fi) e S*dX, s,s'e(0,7r), s.t. 
/, sins' 

o' = tant^ = ^ , r =Acoss, t =— Acoss, 

sm s 

(y',Ai') = Asins'cxp(s'iJi,J(y,/i), (y",/i") ^ -Asinsexp(si7i,J(y,/l)} 

ur+(A)ur_(A), r±(A) = {(a,y,2/,±A,TA,0,0) : a G (0, oo), y e aX} 

The sets T±(A) are, for fixed y, integral curves of both vector fields, and they 
appear separately only because we used the parameterization of Melrose- Zworski. 
The smooth structure near T±(A) follows from the flowout description, but is not 
apparent in this parameterization; we discuss it below while describing the closure 
of L°(A). 

The closure L{X) of L°{X) is ^'^T*X^ is clear from the above description; it is 
(4.8) c\L = L{X)U FxU F-x 

where 
(4.9) 

L(A) ={(0,y',y",r',r",/.',/i") : 3{y,ft) e S*dX, e [0,^], 

(sins)^ + (sins')^ > 0, s.t. 

sins' , ^ / // 

a = t&nO — — , T = Acoss , t = —Acoss, 

sin s 

{y',fJ'') = Asins'exp(s'i7i^)(y,//), (y",/i") = -Asinsexp(si7i,J(y,/i)} 
UT+(A)ur_(A), and 
Fx — {((T, y', y'\ —A, —A, 0, 0) | 3 geodesic of length n connecting y\ y"}. 

Note that the requirement (sins)^ + (sins')^ > just means that s and s' cannot 
take values in {0, tt} at the same time. The set L{X) \ L°{X) comprises those points 
where one of s, s' takes values in {0,7r} while the other lies in (0, tt). The sets 
r±(A) in ( [f.S| ) comprise the limit points where s and s' converge either both to or 
both to TT, whilst F±x comprise the limit points as s — > and s' — > tt or vice versa. 

The smooth structure near T±{X) becomes apparent if we note that near r' — A, 
t" = -A, a e [0, C) where C > 1, L{X) is given by 

{{cJ,y',y",T',T",^i',^i")■. 3{y,fi)eT*dX, |mI < a G [0, C) s.t. 

(4.10) r' = A(l-|a/i|2)l/^ r" = -A(l - H^)!/^, 

{y',fi') = Xexp{f{a^i)Vh){y,a^i), {y",fi") = -XcMfit^)Vh)iy, 

where f{fi) — arcsin(|/i|) is smooth and /(O) = 1. Thus, the differential of the 
map 

(4.11) (y, ^i) ^ -XcMf(p)Vh){y, /i) = (/, Ai") 
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is invertible near /i = 0, so it gives a difFeomorphism near — 0. Hence, a and 
(y", /i") give coordinates on L(A) in this region, so L(X) is smooth here. Away from 
r+(A), coordinates on L(A) are cr, y", /i" and s. 

In the coodinates {y, jl, s, s'), the vector field V; is given by sins'^s' and Vr is 
given by sin s9s- The intersection of L(X) and N* diagjj is given by {s — s'}. Thus 



L+(A) is given by {s < s'}. On L+(A), t = t' + err" by (2.24), so 



sin(s - s') 

T ~ A < on i+(A). 

sin s 

Thus, any distribution in /'"(A^* diagj^, L+(A)) satisfies condition ( |4.2| ). 
We also define 

(4.12) L*{\)^{{e,y',y",-\-\Q,Q): y' ,y" e dX, 9 e /2]}, 
so i'^(A) is a Legendrian submanifold of ^^T*^^X'^, and 

(4.13) c\Lr]L*{±X) ^ F±x. 
Proposition 4.1. T/ie pair 

(4.14) i(A) = (i(A),L#(A)UL#(-A)) 

is a pair of intersecting Legendre manifolds with conic points. 

Proof. We must show that when the set {tq | t > 0, g G L^(A)} is blown up inside 
s<i>rp*j^2^ the closure of L{X) is a smooth manifold with corners which meets the 
front face of the blowup transversally. Let us restrict attention to a neighbourhood 
of L^(A); the case of L'^(— A) is similar. Consider the vector field Vi + Vr. By 



( ^.25| ) and ( ^.26| ), in E(iJ, - A^) n E(iJ,. - A^) this is given by 

2{t' -T")ad,+2T' ^l' ■d^,,+2T" ^l" ■d^.n+d^^h' ■dy,-dy,h' ■d^,,+d,,nh" -dyn^dy.h" ■d^.n 

This is equal to — 2A times the b-normal vector field fi' ■ d^i + fx" ■ d^n plus a sum of 
vector fields which have the form pV , where p vanishes at L*(A) and V is tangent 
to lb and i*(A) (all considerations taking place inside S(-ff/ — A^) n Tj{Hr — A^)). 
Thus, under blowup of {tq | t > 0, g G L*(A)}, Vi + Vr lifts to a vector field of the 
form 

(4.15) Vi + Vr^ -2\sds + sW, 

where W is smooth and tangent to the boundary of L{X), and so dividing by s 
yields a nonvanishing normal vector field plus a smooth tangent vector field. As 
above, such a vector field has a continuation across the boundary to the double of 
L(A) (across the front face) as a smooth nonvanishing vector field. This holds true 
smoothly up to the corner with lb, so L{X) is a smooth manifold with corners. □ 

4.4. Smoothness of symbols. In the next stage of the construction, wc solve 
away the error i?2(A) which is microsupported in the interior of L+(A). This in- 



volves solving the transport equation globally on i+(A). In view of Proposition 4.1 
we can expect the construction to involve Legendrian conic pairs with respect to 
(i+(A), L'^(A)). In order for the symbol to be quantizable to such a conic pair, we 
need to show regularity of the symbol on L+(A), so that it lies in the symbol space 
of the exact sequence from Proposition That is, the symbol of order j — 1/2 
on i+(A) should lie in 

(4.16) p:f-Vr/'"Vf-'^/'"'C-(L+(A);f^,i/^«5t^-i/2l(^+(A))). 
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(We will ignore the symbol bundle S'^-'^^/^l in the rest of this section.) 

To do this, we observe that the symbol on L+{X) automatically satisfies the 
transport equation for the right Hamilton vector field. To see this, let G2(A) be the 
approximation to R{X) constructed so far, with 

and 

'^WF(£;2(A)) C L+{X) n {r < -c} 

for some c > 0. Consider applying Hr — X^ to E2{X). Since Hi and Hr commute, 
and HiKu = HrKi^, we have 

{Hi - X") {{Hr - X^)G2{X) - Id) = {Hr - X^) {{Hi - X^)G2{X) ~ Id) 

- {Hr - X^)E2{X). 

We claim that '*'=WF((i7r - A^)G'2(A) - Ku) is contained in {r < -c}. For if there 
is a point where r > — c, then by (|2.8| ), the maximal bicharacteristic ray in 

(4.17) ^{Hi ~ A2) \ -WF((i?, - X^){{Hr - X^)G2{X) - Km) 

= J:{Hi - A2) \ '-WF{{Hr - X'')E2{X)) 

lies in ^'^WF((i7r ^ A^)G2(A) — i^id)- Such rays always propagate into r > 0. But 

^'=WF(G2(A)) C {r < 0}, ""WF(i^id) C {r < 0}, 

so this is impossible. Consquently, {Hr~X^)G2{X) — KM has no scattering wavefront 
set for {r > — c}, and so the symbols of G2(A) must obey the right transport equa- 
tions in this region. By cutting off the symbols closer and closer to the boundary 
of L+(A), we see that the right transport equations must be satisfied everywhere 
on L+{X). 

Let us examine the form of these transport equations at the boundary of L+(A). 
Near lb, we have coordinates (j/",/i",cr) near T+(A) and {y",p,",a,s) away from 
T+(A), which are valid coordinates for cr < 2, say. The situation near rb is similar 
so the argument will be omitted. 

In either set of coordinates, the left vector field, restricted to L+(A), takes the 
form 

Vi = 2T'ad^. 

Also, by Lemma 2.2, the subprincipal symbol of Hi — X^, which is equal to the 
subprincipal symbol oi H ~ X^ in the singly-primed coordinates, vanishes where ^' 
vanishes, and /i' = at lb on L+(A). Therefore, by (t3.5|), the transport equation 
for the symbol of order —1/2 takes the form 

(4.18) [-^{Cv,-nT')+af)ao\'^dy"d^l"\^ ^0, / G C°°(L+(A)), 
near T+(A), or 

(4.19) ( - i{Cv, - nr') + a/)ao|^yd2/"dA"l^ - 0, / G C°°(i+(A)), 
away from r+(A), which gives an equation for oq of the form 

(4.20) -zT'(a, + /)(a-"/2ao) = 0, / G C°°(i+(A)). 
This shows that (T~"/^ao is smooth across a — Q. 
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To show regularity near L*(A), we use the fact that the symbol satisfies both 
the right and left transport equation. We take the sum of the transport equations 
that obtain when we use the total boundary defining function x' for Hi , and x" for 
Hr- The right transport operator with respect to x" takes the form 

- 71t") +Psub 

However, by (|3.2[) the symbol written in terms of x" is equal to (a;"/x')^^/'^""/^ 



times the symbol written in terms of x' . Since we are writing the sym bol in terms 
of x', we must include a factor cr~^/^~"/^ to be consistent with ( 4.19 ). This gives 
an equation for ao of the form 

(4.21) -nr'O +ftub(y^M^r'0)(a-l/2-"/^ao|^^d2/''dA''|^) =0. 

In view of the term —2t" ad^ in the formula (2.26) for V^, and since psub vanishes 
at s = since /i" = there, we get an equation for oq 

(4.22) (K+T" + s/')ao = /' G C°°(L+(A)). 



Combining with the left transport equation gives an equation which, using (4.15) 
and the fact that r' = t" = — A at s = takes the form 

71 — 1 ~ 

2\{sds — -^sW + sf)ao = 0, 

where W is tangent to the boundary of i+(A) and / is smooth on i+(A). This 
may be written 

(a. + W^ + /)(.s-("-i)/^ao)=0, 
This together with ( |4.20| ) shows that ao G (j'"/'^s^'"-^'^/'^C°°{L+{X)). 

It follows that there is a Legendre distribution /~^'^~'^~'~2~ (L_|_(A), L+(A)) 
which has the correct symbol of order —1/2 at L+(A). Thus it solves the equation 

(4.23) {Hi - X')Q - E{X) e l'-'^-'^^^{XlL+{\),L*{Xy,^''ni), 

where (ri + 3)/2 is the order of vanishing at lb and (ri — 1)/2 is the order of vanishing 
at rb. The order of improvement at lb is two since not only is the Legendrian Gi of 
Definition 2.5 at lb characteristic for Hi — X^, but the transport operator for symbols 
of order (n — l)/2 vanishes, so we automatically get two orders of improvement 
there. At rb however we can expect no improvement. As shown above, Q will 
automatically satisfy the equation 

[Hr - X')Q - E2{X) e /l.f ;^'^(X2, L+(A), L#(A); ^*r!^). 

Let us assume by induction that we have a kernel which solves the left equation 
above up to an error in 

(4.24) {Xl,L+{X),L*{X)f'^n^) 
and hence the right equation up to an error in 

(4.25) {Xl,L+{X),L*{X)f^n^) 

We wish to improve this by one order at L+(A). To do this, we choose Qk E 
/''■~5'T-;-T-.T-(L_|_(A),L+(A)) to have the symbol of order k — 1/2 on L+{X) 
which solves the left transport equation (and therefore the right transport equation) 
on L+{X). We need to investigate the regularity of this symbol to see if it extends 
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to a Legendrian conic pair. The argument is very analogous to the one above, but 
now we have error terms on the right hand side. In the first region, after removing 
the half-density factor, we get an equation of the form 

(4.26) -2iT'{ad„ + i-- + k)+af)ak^bk, / G C°°(L+(A)). 

The term bk comes from the error to be solved away. Since the error term is of 
order fc + 1/2 at L+(A) and order {n + 3)/2 at lb, bk € cr"/2-'=+^C°°(L+(A)). This 
shows that Ofc G cr"/^^'^C°°(i+(A)), as desired. Similarly, in the second region, near 
the corner Ibni^(A), by combining the vector fields Vi +Vr we get an equation of 
the form 

(4.27) ~2iT"(^-sds + sW+{--^ + k) + saf)ak^bk, feC^{L+{X)) 

with bk again the error to be solved away. To calculate its order of vanishing at 
s = 0, consider the transport equation for symbols of order (n — 2)/2 at L'^{\). 
Noting that the subprincipal symbols vanish identically on _L^(A), the left transport 
operator is 

'~i{Cvi - t') 

whilst the right transport operator with respect to x" is 



To write this with respect to x' we must conjugate by cr (by (3.2)). In view of the 
term — 2r"(T9cr, this changes the operator to 

-^{Cv,.+r"). 

The sum of these two operators vanishes on L*(A) so actually the right hand side 
in (4.27) comes from a term in {X^, L^{X), L^{\);^'^fl^). From 

Proposition 3.5 we see that bk e s("^^)/^^'^+V"/^^'^C°°(i+(A)), one power in s 



better than might be expected. This shows that Ofc € s("-i)/2-'=cr"/2-'=C°°(L+(A)), 
as desired. Therefore, one can find a Legendre conic pair with symbol of order 
A: — 1/2 on L-f (A) equal to Ofc which solves away the error term of order k + 1/2 
at i+(A). This completes the inductive step. By asymptotically summing these 
correction terms, we end up with an approximation Gs{X) to the resolvent kernel 
with an error £;3(A) in J^;'^''^(X2,L#(A);^*r!5). That is, we have solved 
away the scattering wavefront set of the error term at L^{X) completely. 

4.5. Solving away outgoing error. In the last step of the construction of the 
parametrix, we solve away the error to infinite order at bf and lb. We begin by 
considering the expansion at rb. By construction, the parametrix G3(A) constructed 
so far has an expansion at rb 

iA/a:"^//(n-l)/2^^,/J^ 
j>0 

where gj{z',y") ■ v' e /-"/4-J>/'*-i/2-j(G'j^„(a), G«(A))(X;^'=f2^), Gy"(A) is the 



G3(A) ~ e'^/" x"^""'^'')_Ja;"^5j(z',y") • v' ■ u" 

j>0 



fibre Legendrian of Definition 2.5 and y" is regarded as a smooth parameter. The 
factors i^' , v" are the Riemannian half-density factors on X lifted to X"^ via the 
left and right projections, respectively. We will ignore the half-density factors from 
here on; since A{a ■ v' ■ v') = A{a) ■ v' ■ v' , this only has the effect of changing 
= Z\ + P to Z\ + P' for some P' with the same properties as P. 
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The error term after applying H — to G3(A) has the form 

(4.28) EsiX) ~ e'^/""a;"^""''/'^a;"^e,(z',y"), 

j>0 

where Cj S ^''■"'''^''^^"■'e''*'/^ C°°{X x Y). Again we regard y" as a parameter. Thus 
we have 

(4.29) {H - X')g, = e, 

Consider the problem of solving away errors of the form ej , to infinite order at 
bf (of course we cannot solve the errors away exactly without begin able to solve 
{H — X^)u = f exactly, which we cannot do until we have constructed the resolvent 
kernel!). If we apply — to a series of the form 

(4.30) e*^/^x("-i)/2-fc^2.j^^.^ 6jeC°°(X), 

j>0 

we get a series of the form 

(4.31) e*^/^x("+i)/2-'=^2;^'cj, CjGC°°(X), cq ^ 2iXkbo. 

3>0 



Thus, we can add to gj a series of the form (4.30) to solve away the powers greater 



than (n + l)/2, but the power (n + l)/2 presents a problem (without introducing 



logarithmic terms), because of the vanishing of k in (4.31), unless the coefficient of 
the term happens to be zero. We need the following results. 

Lemma 4.2. If g £ e'^/'=x'^''-^^/^-''C°°{X), k^l,2,..., satisfies 

(iJ-A2)g = x("+i)/2e*^/^C°°, 

then 

(H-A2)5 = a;("+3)/2g«A/a;^oo^ 



Proof. It follows inductively using (4.31) that the coefficient of order (n + l)/2 — ? 



vanishes, for / = fc, A: - 1, . . . . Thus, actually g G e'^/'=x'-'^-'^^^^C°°{X). Then ( p3l|) 



shows that the next coefficient also vanishes. □ 

Corollary 4.3. The same result holds if the condition g e g*>/a:2;("-i)/2-fe(joo^j^^ 
is replaced by g £ /P.("-i)/2-"/4-fcj-^^ ^^iJ-j ^ j^^ggjidj-g conic pair 

Proof. Apply the above argument to the symbol at G'*. □ 



Thus, for each j, we can modify gj by a series of the form (4.3C) until the error 
term is of the form x'''"+^'^^^e'^/^'C^{X x Y). Then applying the Corollary to g^, 

we find that unsolvable term of order x'^^^^^^'^ vanishes. Therefore, we can solve 
away the ej to infinite order at bf. Thus, we may assume that our error in i?3(A) 
vanishes to infinite order at the corner bf n rb. 

Next, we solve the error away at L"^ . This involves solving the transport equation 



(4.32) «A(a9. + (i+j))aj = 

The equation for aq then is 



by 
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and bo is rapidly decreasing at rb and is in a'^/^C°°{L'^{X)) at lb. There is a unique 
solution which is rapidly decreasing at rb and in cr~^/^C°° at lb. We can thus find 
a correction term which reduces the error to /("+2)/2;("-+3)/2,(n-i)/2^2,#(A)), with 
infinite order vanishing at bf n rb. Inductively, assume that we have reduced the 
error to /»/2+fe;(n+3)/2,(n-i)/2(^#(;)^))^ .^yj^h infinite order vanishing at bf n rb. The 
transport equation for is then 

(crdcr + ^ + fc^flfc = bk, 

where inductively, bk is rapidly decreasing at rb and is in (7^/^~*^C°°(L*(A)) at lb. 
There is a unique solution rapidly decreasing at rb and in a^^/^^''C°°{L'^{\)) at lb. 
A Legendre distribution in /»/2+fc-i,(n-i)/2,(«-i)/2(-2,#(A)) with a-j as symbol then 
reduces the error to /"/2+'''+iK"+3)/2,(n-i)/2^^#^^-j^^ with infinite order vanishing 
at bf n rb, so this completes the inductive step. Taking an asymptotic sum of such 
correction terms yields a parametrix 6*4 (A) — G(A) leaving an error which is the 
sum of a term supported away from rb of the form 

with a smooth and rapidly decreasing at bf, plus a term supported away from lb of 
the form 

at rb, with b rapidly decreasing at bf. The error at lb can be solved away using 
(pel)- dial]), leaving an error term -B4(A) which can be expressed on the blown- 



down space as 

S+(A) = e'^/""x"^"~'^/'5(z',z") 

with b smooth on and rapidly decreasing at x' — 0. Such an error term is 
compact on the weighted l? space x^lP'iX') for any ^ > 1/2 (where I? is taken with 
respect to the metric density) . This completes the construction of the parametrix. 

5. Long range case 

The case of long range metrics and long range perturbations, P € xDiSl^{X), 
requires only minor modifications in the parametrix construction until the last step 
(removing the outgoing error). In particular, there is no change in the construction 
of the pseudodifferential approximation. In the intersecting Legendrian construc- 
tion, as well as solving the transport equations on L°{X), the only difference is in 
the structure of the subprincipal symbol, which no longer obeys Lemma ETsI Thus, 



the arguments in section 4.4 and section 4.5 have to be modified. Let q denote the 



boundary subprincipal symbol of H . Notice that in the gravitational long range 
case, g is a constant, but in the general long range case, q is an arbitrary smooth 
function of y which is a quadratic on each fibre of K over Y. Let qi and qr denote 
the lift of q to X^ via the left, respectively right, projection. 



Let us now discuss the necessary modifications to sections 4.4 and 4.5. Equation 
( [4.1S| ) becomes 

(5.1) -iijO-Vi - nr' + iqi)ao\^dy"d^i"\^ = 0, 
near T+(A). Thus, ( [4.20| ) is replaced by 

(5.2) -zT'(9, + /)(a-"/2-lta„) = o, / e C°°(L+(A)). 
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Thus, now we conclude that in this region ao is of the form 

(Ti+'^C°^{L+{X)). 

Next, in the second region, at the corner lbnL'*(A), the right transport equation 
( [4.22| ) becomes 

(5.3) {Vr+r" + iqr + sf)aa = /' G C°°(L+(A)). 
Adding this to the left transport equation yields 

^ „ n — I Ql Qr ~^ 

2X{sds - + + = 0, 

which now gives that oq is of the form 

at+*S:s'^+'^C°°(L+(A)). 

Combining this with the other similar results at rb and the interior of (A) , we 
deduce that 

In the general long range case, the dependence of q on y, and its appearance in 
the exponent of the boundary defining functions pib, etc., means that differential 
operators from the left factor, acting on a Legendre function with principal symbol 
ao, introduce logarithmic terms. For example, in a neighborhood of lb in L+(A) 
the error term bk in (4.26) for fc = 1 will take the form 

6i ^a^+'^((loga)2c2 + logaci+co), c, G C°°(i+(A)), j = 0,l,2. 
Then the transport equation for ai takes the form 

(5.4) (^-t{Cv,+i-n + 2)T' + iqi)+afyi\^dy"dp"\^ =0, / G C°°(L+(A)), 

n 

(5.5) -ir'{ad,^- + l-iqi/i2\) + af)a,=bu / G C°°(L+(A)). 
Hence, near lb but away from L'^(X), ai will take the form 

ai =at-i+'lt ((log a)24 + log a c'l + c^), c;. G C°°(i:+(A)). 



A similar discussion works at the other boundary faces of L+(A), with up to 
quadratic factors in each of logpib, logprb, log^^, and can be repeated (with pro- 
gressively higher powers of logarithms) for all 's. 

Since the most important long-range case is the gravitational case where the 
subprincipal symbol is constant, and since it makes the discussion more transparent, 
in what follows we make the assumption that 

A and P are of long range gravitational type, 

which implies that q is constant. Let 

q 

"^2A- 

The point is that in this case the powers of pib, etc., above are constant, thus no 
logarithmic factors arise when we apply _ff — A^ to such Legendre functions. Then 

ak e pr/'+*>r/'+'""Vf "'^^'+''""'C°°(^+(A);f^^/' 5[-i/2l), 
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and asymptotic summation gives an outgoing error 

^+(A) e e'V-'e'V-"2,"(»-i)/2+-2.'(«+i)/2+-coo(^2.sc^i)^ 



Since a is a constant, ( 4.3C )-( 4.31 ) are still true (fc need not be an integer; it suffices 



that it is a constant), except that now 

Co = i2iXk + 2Xa)bo. 

Since now we are taking k = I + ia, where / is an integer, we can solve away the 
series, provided that the coefficient of the / = term vanishes, which is assured just 



as in Lemma 4.2. The rest of the argument requires only similar modifications as 
compared to the short-range case, so we conclude, as there, that we can modify the 
parametrix to obtain an error term of the form 

with h smooth on and rapidly decreasing at x' = 0. 

6. Resolvent from parametrix 

In the previous two sections, we constructed a parametrix G{\) for i?(A) which 
satisfies 

{Hi-\^)G{\)^Ku + E{\), 

where E{\) has a kernel which is of the form e^^/^"a;"^""^^/^a;'°°C°°(X^; '^'^f^s). 
Thus, it is a Hilbert-Schmidt kernel on x^lJ^{X) for every / > 1/2, and in particular 
is compact. In fact, we also see directly from its form that E{X) : x^L?'{X) 
C°°{X) for / > 1/2. Crude estimates (such as Schur's Lemma) show that G(A) acts 
as a bounded operator from x'L^ to x~''L^ for large enough / > 1/2; more refined 
estimates, which we do not need here, show that in fact this is true for any I > 1/2. 
Thus, the equation above becomes an operator equation 

{H~X^)G{X) =ld+E{X) 

from .t'L^ to x^^L^. 



6.1. Finite rank perturbation. To correct G{X) to the actual i?(A), we must 
solve away the error term E{X). Thus, we would like Id+£'(A) to be invertible. 
However, this is certainly not necessarily the case as things stand; if for example we 
modified G(A) by subtracting from it the rank one operator G(A)(0)(0, •), for some 
(j) e {X , then the modified G(A) would be microlocally indistinguishable 
from the old one, but would annihilate 0, so the modified Id +i?(A) would not be 
invertible. 

Since Id+_E(A) is compact, it has a null space and cokernel of the same finite 
dimension N . To make Id +E{\) invertible, we try to correct G(A) by adding to it 
a finite rank term 

N 

(6.1) 5^</..(x2V^,-)- 

i=l 

Here (, ) denotes pairing, 7/;^ should lie in L^, and the factor of x^^ is included 
to ensure that it acts on x'L^. We require that 0^ are in x^'i^ so that ( |6.lD is 
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bounded from x'L^ x^'L^. We wish to choose and i/^i so that 

N 

ld+E{\) + Y,{{H - A2)0,)(x2V., •) 



is invertible. This is possible if we can choose x'"ipi to span the null space of 
Id+i?(A) and {H — \^)(f)i to span a subspace supplementary to the range. Note 
that if (Id +E{X))u — and u S a;'L^, then u — —E{X)u, so the mapping properties 
of E{X) imply that u £ C°°{X). Thus, we automatically have tpi G C^{X) above. 
To proceed, we need the following lemma. 

Lemma 6.1. Let I > 1/2. Then the image of H — }? on the sum ofC°°{X) and 
the range of G{X) applied to C°"{X) is dense in x^V^. 

Remark 6.2. Note that {H-\^)G{X)g ^ {ld+E{X))g e C°°(X) if g G C°°(X), and 
for u e C°°{X), {H - \^)u e C°°(X), so the image of F - on the space in the 
statement of the lemma is a subspace of C°°(X). 

Proof. To proceed, we give the proof for short range H; the proof for long range H 
requires only minor modifications. 

Let be the subspace of x^L^ given by the image of if — A^ on the sum of 
C°°{X) and the range of G{\) applied to C°°{X). If M. is not dense, then there is a 
function / G x^ L"^ orthogonal to M. Since for u G C°°{X) imphes {H — \^)u G M, 
f satisfies 

{x-^f,x-\H - X^)u) = Q VmgC°"(X) 

(6.2) ^{{H-X^)x-^^f,u) = {] VugC°°(X) 

^{H-X^)h^Q, h^x-^^f. 

where we used that H is symmetric on C^{X). On the other hand, G(A) maps 
C°°(X) ^ a;("-i)/2e*-^/^C°°(X), andforanyg G C°^{X), {H-X^)G{X)g = {ld+E{Xj)g G 
C°°{X), hence {Id+E{X))g G M. In addition, E{X)*, with kernel E{X)*{z',z") = 
E{X,z",z'), maps C-°°{X) x^'^-^^l'^e-'^I^C°°{X), so we have 

{x-'f,x-\H-X^)G{X)g)^{x-'f,x-\\A^E{X))g)^^ VgGC°°(X) 

(6.3) ^((Id+£;(A)*).T-27,g)^0 VgGC°°(X) 

^ (Id+£;(A)*)/i = 0. 

\ih = -E{X)*h, then h has the form x''"-^^/'^e-'^/^C°°{X), ie it is incoming. A 
standard argument then implies that h = 0. Let h — x'-'^^^^^^e^^^^^hoiy) + h, 
where h G a:("+i)/2e~'^/"^C°°(X). Green's formula yields 

= / h{H - X^)h- h{H - X^)h 
J X 

(6.4) r _ _ r 

= 2a lim / hx^da:h-hx^d.^h^2iX \ho{y)\^, 

{x=e} Y 



so Hq = 0. It then follows iteratively from ( 4.30 ) and ( 4.31 ) that the expansion of 
h at the boundary of X vanishes identically, that is, that h G C°^{X). Finally a 
unique continuation theorem, see e.g. Chapter XVII], shows h = identically. 
This means that A4 is indeed dense inx^L^. □ 
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Thus, we can choose the 0^ G x("-i)/2e*-^/^C°°(X) above so that {H ~ e 
C°° (X) span a supplementary subspace of range Id +i?(A) . The modified parametrix 
then satisfies 

[H - \^)Gr-,{\) =ld+E^{\), 

where £'5(A) has the same form as E{X) but in addition is invertible on x^Lp' for aU 
I > 1/2. 

6.2. Resolvent. The inverse Id+5'(A) of Id+£;5(A) is Hilbert-Schmidt on x'L"^ 
since this is true of £'5(A). Moreover, since 

S{X) = ~E5{X) + E^iXf + E5{X)S{X)E5{X), 

it is easy to see that also S{X) e e*^/^"a;"'""^^/\'°°C°°(X^; ^^^=17^). Our solution 
for the kernel i?(A) is then 

R{X) = G5(A)(Id+5(A)). 
It is not hard to show that G5(A)5(A) has the form 

^,(n-l)/2^„(„-l)/2^a/x'g»A/x"^oo(;^2.sc^i)^ 



so i?(A) has the desired microlocal regularity (4.4). 

Remark. Lemma 6.1 directly shows the absence of positive eigenvalues of 
H. Suppose that {H = and u e x'H'^iX) for some s > -1/2. This 

would certainly be the case of an eigenfunction since H has elliptic interior symbol, 
so u would lie in H^{X) for every k. We need to show that for all functions 
g e a;("~i)/2e^-^/"^C°°(X) the equation 

(6.5) {[H~X'')u,g)^{u,[H-X'')g)^Q 

holds. Indeed, this implies that u is L^-orthogonal to the image of _ff — A'^ acting on 
^{n-i)/2^i\/xQoo(^-^^^^ or equivalently that x^'-u e x^L'^ is orthogonal in x'L^ to the 
image oi H — X^ acting on x^"'~^'>/'^e'^^/^C°^{X) in x^L'^. Then Lemma 3J shows 
that u = 0. 

To deduce (U) for g e x^'^-^'^/'^e^^/''C°^{X), let e C°=(M), identically 1 on 
[l,oo), near the origin. Then 

{{H - X^)u, g) = Imi ((i/ - A^)^, cl>{x/t)g) = Ihn (m, {H - X^)ct>{x/t)g) 

(6.6) = Imi ^{x/t){H - X^)g) + Inn (u, [H, ^{x/t)]g) 

= {u, {H - X^)g) + lim (u, [H, ^{x/t)]g). 

Note that [H,(f){x/t)] is uniformly bounded (i.e. with bounds independent of t) as 
a map x'i?^ x^^^L^, and in fact [H^4){x/t)] — > strongly as i ^ 0. Applying 
this with I — —1/2 — e , e > sufficiently small, we see that the last term goes to 
as t ^ 0, proving (6.5). 

6.3. Analytic continuation. It is not hard to show that the kernel R{X) con- 
structed above continues analytically (as a distribution on X"^^) into Im A > 0, Re A > 
0. We complete the proof of Theorem |l . l| by showing that this analytic continuation 
coincides with the outgoing resolvent, R{X^), for ImA > 0. 

Our parametrix is defined as an asymptotic sum of symbols, which is really a sum 
with cutoff functions inserted (see Proposition 18.1.3] for an explicit construc- 
tion). The cutoffs depend on C'^ norms of a finite number of symbols and ensure 



32 



ANDREW HASSELL AND ANDRAS VASY 



that the sum converges in C*^ for all k. If the symbols are holomorphic in A then 
the C*"' norms may be taken uniform on compact subsets of A. Since holomorphy 
is preserved under uniform limits, we need only check that the phase and symbols 
analytically continue in some explicit parametrization of the Legendrians. 

It is standard that the pseudodifferential approximation Gi (A) analytically con- 
tinues. Blowing down sf, we solve away the error as an intersecting Legendrian, 



see Section 4.2. Let cj) he a, local parametrization of the Legendrian L{1) near 
L(l) n N* diag^j. Then it is easily checked that the phase function 

k<j) + s(A - k) 

locally parametrizes (A''* diagj,, L"'"(A)). Since the variable s takes nonnegative val- 
ues, the function gtii':<P+si^-k))/x continues to ImA > 0. 

Away from sf, the value of r is strictly negative on the Legendrian, and so the 
phase is of the form 

where cf) is positive on the Legendrian, independent of A. By restricting the support 
of the symbol sufficiently, therefore, we may assume that (f) is positive everywhere 
on the support of the integral. Thus this also analytically continues to the upper 
half plane with uniform bounds. 

The symbols are defined by iteratively solving transport equations of the form 

where 6o = 0. These equations are solved along L°(A), i.e. if we consider am- 
plitudes in an explicit parameterization of the Legendrian, then along the critical 
submanifold 6*0 = {(0,y,u) : du4> — 0}, where X(j)/x is the phase function as 
above. Note that is independent of A, and it is identified with L°{X) via the 
map Cj, 3 (0, y, u) i— > {0,y,d(^x.y){X4'/x)\{o,y,u))- Along the transport equation 
becomes an ODE whose coefficients depend on A polynomially, since the only A de- 
pendence of the coefficients arises from this identification map, and Hp, Psub are 
polynomial in the fiber variables. Thus, the solution aj of the transport equation, 
as a function on D x C^, D a neighborhood of the positive real axis in C = C\, is 
holomorphic (in A), provided bj is (here we identify with L°(A)). Note that the 
bj's arise because solving the transport equations only guarantees that the 'error 
term' E^i^X), arising from the application oi H — X^ to i?3(A), is one order lower 
than expected, so for each A, 6j|c^ depends on a^, i < j near C^, and not just on 
ai\c^- (In fact, bj\c^ depends on a finite number of terms of the Taylor series of 
Ui, i < j, at Ctf,.) To ensure that the bj are holomorphic in A, we define the near 
D X C^, rather than at D x C^, e.g. by introducing a local product decomposition 

X [/, t/ c K'^, of the parameter space near C^, and pull-back the a^, first defined 
on Z? X C^, by the projection. Then, having constructed ai, i < j, bj will be holo- 
morphic in A near, hence at, I? x C^, so aj is also holomorphic at D x C^, hence 
it extends to be holomorphic near D x C^. If we express the amplitudes aj with 
respect to a different parameterization of L°{X), which is still of the form X(j)/x, 
then the new amplitudes dj will still be holomorphic functions of A, so holomorphy 
is preserved at the overlap of parameterizations of different parts of L°{X). This 
completes the inductive argument. 

Therefore, our parametrix constructed above may be assumed holomorphic in 
some set B{e, Aq) n {ImA > 0}, for some Aq > 0. It is easy to see that for non-real 
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A, the parametrix is in the small calculus, since the positivity of implies that the 
exponent of e^^'i'l^ has negative real part, and is therefore rapidly decreasing at bf, 
lb and rb. The finite rank correction may be taken independent of A if we chose 
e > sufficiently small. Then, we have 

(Z\- A2)G'5(A) = ld+£;5(A), ImA>0, ReA>0, 

where all terms are holomorphic in some small open set as above, E'^(X) is invertible 
on x^L^ for all I > 1/2, and off the real axis, G5(A) and Er,{X) are in the small cal- 
culus. Define ld+S'(A) to be the inverse of Id -|-i?5(A) on x'-L^ for some fixed I. By 
the symbolic functional calculus j^] , for Im A > 0, 5' is a family of scattering pseudo- 
differential operators which is clearly holomorphic. Then i?(A) = G'5(A)(Id +5(A)) 
satisfies {H — A^)i?(A) = Id on cc'L^. But by self-adjointness, and the symbolic 
functional calculus, for ImA > 0, {H — A^) has a pseudodifferential inverse on L^. 
Since -R(A) is a bounded operator on for Im A > it must be the inverse. There- 
fore we have shown the inverse on the real axis constructed above continues as a 
Schwartz kernel to the upper half plane and agrees with the resolvent there. This 
completes the proof of Theorem 1.1. 

Remark 6.3. The only place where we used that Im A > is to make our parametrix 
act on, and its error compact on, weighted Sobolev spaces. Namely, in the last step 
of the construction, i.e. when we add a finite rank perturbation to remove the error 
E{X), we need E{X) to be a compact operator on x'i^ for / > 1/2. However, the 
kernel of i;(A) is of the form e*^/^"a;"'""^^/\'°°C°°(X^; "'=17^), and for ImA < the 
real part of the exponent is positive, so the kernel of E{X) is not even a tempered 
distribution on X^. In particular, it does not even map C°°{X) to C~°°{X). The 
same statement holds for G(A) as well. 
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